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INTRODUCTION 


Research continued on Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell. 
The motivation for this project is the planned utilization of advanced composite materials in the 
fuselage for large transport aircraft. In particular, the focus of this activity is the effect of cabin 
pressurization on the stiffener-to-skin joint. The design of stiffener-to-skin joints is one of the 
major technology issues in utilizing graphite/epoxy composites in the fuselage of a large transport 
aircraft. The manner in which the loads are transferred in the stiffener-to-skin joints under inter- 
nal pressurization is important for determining the load capacity of these joints. 

The objective of this project is to develop analyses of an orthogonally stiffened composite 
cylindrical shell subjected to internal pressure. These analyses are used to study the distribution 
of the interacting loads between the shell and stiffeners, and to study the pillowing of the shell, for 
a geometry and pressure typical of a large transport aircraft. Primarily the aim is to understand 
the fundamental mechanics of the load transfer in the vicinity of the shell-ring-stringer joint. Sec- 
ondly, these analyses can be used in parametric studies of joint response, and perhaps for design. 
A potential benefit of such an analysis/design capability is to use fewer expensive fasteners in the 
graphite/epoxy fuselage. Where fasteners are required in a graphite/epoxy structure, aluminium 
fasteners cannot be used because of galvanic corrosion to the metal. More expensive fasteners, 
like titanium, are required to avoid corrosion. 


RESEARCH ACCOMPLISHED 

A conference paper was presented, and an extended abstract of it appears in the conference 
proceedings. The citation for the abstract and presentation are given below. 

• Johnson, E.R., and Rastogi, N., “Influence of an Asymmetric Ring on the Modeling of an 
Orthogonally Stiffened Cylindrical Shell,” Proceedings of International Conference on 
Composites Engineering ICCE/1, David Hui, Editor, International Community for Com- 
posites Engineering, August 1994, pp. 237 & 238. 

• Johnson, E.R. (speaker), and Rastogi, N., “Influence of an Asymmetric Ring on the Mod- 
eling of an Orthogonally Stiffened Cylindrical Shell,” International Conference on Com- 
posites Engineering, Sheraton Hotel Downtown, New Orleans, Louisiana, August 28-31, 
1994, Session 13f: Composite Structures 1. 

Work continued on the effect of a ring, or frame, with an asymmetrical open cross section on 
the response. Since the ring has an asymmetrical section, it twists and bends out-of-plane under 
the internal pressure load in addition to bending in its plane and stretching along its circumfer- 
ence. For the structural repeating unit shown in the Figure below, the asymmetrical section ring 
results in loss of symmetry of the deformation about the 0-axis. (Symmetry about the x-axis is 
preserved since the stringer cross section is assumed to be symmetric.) 

The major new feature incorporated into the analysis was the warping deformation of the 
ring’s cross section due to torsion. This warping deformation is in addition to a previous exten- 
sion of the model to include deformations due to transverse shear in the stiffeners and in the shell. 


Grant NAG-1 -5 37 


2 



z 



Figure: Structural repeating unit of an orthogonally stiffened cylindrical 
shell subjected to a internal pressure p. 


Inclusion of warping was found to significantly change the torsion and out-of-plane bending 
response of the ring, and to change the distributions and magnitudes of the interacting line loads 
between the stiffeners and the shell. For example, the direction of the rotation about the circum- 
ferential axis of each of the structural elements at the shell-ring-stringer joint was change by the 
inclusion of this warping deformation. Also, the circumferential component of the moment 
resultant at the joint due to these interacting load intensities was increased by the inclusion of 
warping deformation. These results are detailed in the Appendix of this report. 
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APPENDIX 


Analysis of an Internally Pressurized Orthogonally Stiffened Cylindrical Shell 

with an Asymmetrical Section Ring 



ANALYSIS OF AN INTERNALLY PRESSURIZED ORTHOGONALLY STIFFENED 
CYLINDRICAL SHELL WITH AN ASYMMETRICAL SECTION RING 


Naveen Rastogi* and Eric R. Johnson^ 

Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24061 

ABSTRACT 

The linear elastic response is determined for an internally pressurized, long 
circular cylindrical shell stiffened on the inside by a regular arrangement of identical 
stringers and identical rings. Periodicity of this configuration permits the analysis of 
a portion of the shell wall centered over a generic stringer- ring joint; i.e., a unit cell 
model. The stiffeners are modeled as discrete beams, and the stringer is assumed to 
have a symmetrical cross section and the ring an asymmetrical section. Asymmetry 
causes out-of-plane bending and torsion of the ring. Displacements are assumed as 
truncated double Fourier series plus simple terms in the axial coordinate to account 
for the closed end pressure vessel effect (a non-periodic effect). The interacting line 
loads between the stiffeners and the inside shell wall are Lagrange multipliers in the 
formulation, and they are also assumed as truncated Fourier series. Displacement 
continuity constraints between the stiffeners and shell along the contact lines are 
satisfied point-wise. Equilibrium is imposed by the principle of virtual work. A 
composite material crown panel from the fuselage of a large transport aircraft is the 
numerical example. The distributions of the interacting line loads, and the out-of- 
plane bending moment and torque in the ring, are strongly dependent on modeling 
the deformations due to transverse shear and cross-sectional warping of the ring in 
torsion. 
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INTRODUCTION 


The design of stiffener-to-skin joints is one of the major technology issues in 
utilizing graphite-epoxy composites in the fuselage of a large transport aircraft 
(Jackson, et ah, 1984). Stiffeners can be attached to the skin by either fasten- 
ers, co-curing, adhesive bonding, or some combination of these methods. Where 
fasteners are required in a graphite-epoxy structure, aluminium fasteners cannot 
be used because of galvanic corrosion to the metal. More expensive fasteners, like 
titanium, are required to avoid corrosion. Hence to reduce manufacturing costs, 
mechanical fasteners can be eliminated in favor of bonded joints. As an example, a 
graphite-epoxy crown panel for the fuselage of a large transport aircraft was recently 
fabricated without fasteners by co-curing the stringers and co-bonding the rings, or 
frames, to the skin (Ilcewicz, et ah, 1992; Swanson, et ah, 1992). Also, the curved 
graphite-epoxy fuselage frames were manufactured by resin transfer molding into 
two-dimensional braided preforms of net structural shape (Jackson, 1994). Clearly, 
the strength of the bond line is a critical issue for these primary fuselage structures 
made from advanced composite materials. The purpose of this paper is to ana- 
lyze the load transfer in bonded stiffener-to-skin joints under cabin pressurization. 
Internal pressure is an important load to consider because it tends to cause peel 
stresses in the bond line which are particularly debilitating in adhesive joints. 

An idealized structural model of the fuselage is analyzed. This configuration 
is a long circular cylindrical shell stiffened on the inside by a regular arrangement 
of identical stringers and identical rings (frames). Periodicity of this configuration 
permits the analysis of a portion of the shell wall centered over a generic stringer- 
ring joint; i.e., deformation of a structural unit cell determines the deformation of 
the entire shell. The stringer is assumed to have a symmetrical cross section , and 
the frame is assumed to have an asymmetrical open section. Asymmetrical open 
section frames are commonly used as transverse stiffeners in the fuselage structure. 
The stiffeners are modeled as discrete beams perfectly bonded to the inside shell 
wall, so that the interacting loads between the stiffeners and shell wall are line 
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load intensities. These line load intensities represent resultants of the tractions 
integrated across the width of the attachment flanges of the stiffeners. 

Mathematical formulations for the linear elastic response presented in this paper 
include the effect of transverse shear deformations and the effect of warping of 
the ring’s cross section due to torsion. These effects are important when the ring 
has an asymmetrical cross section, because the loss of symmetry in the problem 
results in torsion of the ring, as well as out-of-plane bending, and a concomitant 
rotation of the joint at the stiffener intersection about the circumferential axis. 
This stringer-ring-shell joint is modeled in an idealized manner; the stiffeners are 
mathematically permitted to pass through one another without contact, but do 
interact indirectly through their mutual contact with the shell at the joint. Restraint 
of cross-sectional warping, as occurs here in the ring due to contact with the shell, 
is an important contributor to the normal stresses in thin-walled open section bars, 
as was demonstrated by Hoff (1945). Based on transverse shear deformation and 
cross-sectional warping of the ring, four structural models can be defined. The 
simplest model uses non-transverse-shear-deformable theory, or classical theory, and 
neglects warping due to torsion. The most complex model includes both effects. 
Intermediate complexity models occurs for inclusion of one effect without the other. 

For symmetric section stiffeners, the response of the unit cell is symmetric about 
the stringer axis and the ring axis, and there is no rotation of stringer- ring-shell joint. 
Results have been published for the linear response (Wang and Hsu, 1985) and for 
the geometrically nonlinear response (Johnson and Rastogi, 1994) of a symmetric 
configuration subjected to internal pressure. 

MATHEMATICAL MODEL 

An idealized mathematical model is assumed for the semi-monocoque fuselage 
to study the generic characteristics of the response in the vicinity of the stiffeners’ 
intersection. The model is of a very long circular cylindrical shell internally stiffened 
by identical stringers equally spaced around the circumference, and identical frames 
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or rings, equally spaced along the length. In general, the spacing of the stringers is 
not the same as that of the rings. The structure is periodic both longitudinally and 
circumferentially, and the loading is spatially uniform. Consequently, a structural 
repeating unit (or unit cell) can be defined whose deformation determines the defor- 
mation of the entire structure. A typical repeating unit consists of a portion of the 
shell wall centered over portions of stringer and ring as shown in Fig. 1. The radius 
of the middle surface of the undeformed cylindrical shell is denoted by R , and the 
thickness of the shell is denoted by t. Axial coordinate x and the circumferential 
angle 8 are lines of curvature on the middle surface, and the thickness coordinate 
is denoted by z, with —t/2<z<t/2. The origin of the surface coordinates is 
centered over the stiffeners intersection so that — / < x < / and —0 < 8 < 0, where 
21 is the axial length, and 2 RQ is the circumferential arc length of the repeating 
unit. 

The stiffeners are mathematically modeled as one-dimensional elements, or dis- 
crete beams, so that the actions transmitted by the stiffeners to the inside of the 
shell wall are represented by distributed line load intensities. In this paper it is as- 
sumed that the stringer is symmetric about the x-z plane through its centroidal axis 
and the ring is asymmetric. On the basis of the symmetry about the z-axis for the 
unit, only the interacting line load components tangent and normal to the stringer 
are included in the analyses. The shell-stringer interacting force components per 
unit length along the contact lines are denoted by A X3 (x) for the component tangent 
to the stringer, and A 25 (x) for the component normal to the stringer. However, due 
to an asymmetrical cross section ring, the components of line loads between shell 
and the ring consist of three force intensities and two moment intensities. The three 
shell-ring interacting force components per unit length along the contact lines are 
denoted by A rr (0) for the component acting in the axial direction, A e r {8) for the 
component tangent to the ring, and A 2r (0) for the component normal to the ring. 
The two shell-ring interacting moment components per unit length along the con- 
tact lines are denoted by Ao r (8) for the component tangent to the ring and A zr (8) 
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for the component normal to the ring. These interacting loads acting in a positive 
sense on the inside surface of the shell are shown in Fig. 2. The purpose of the 
analysis is to determine these distributed line load intensities and also, to examine 
the differences in their magnitudes and distributions for the four structural models 
described earlier. 

For all the structural models, the linear elastic response of the repeating unit to 
internal pressure is obtained by utilizing Ritz method and the principle of virtual 
work applied separately to the shell, stringer, and ring. The virtual work functionals 
are augmented by Lagrange multipliers to enforce kinematic constraints between the 
structural components of the repeating unit. The Lagrange multipliers represent 
the interacting line loads between the stiffeners and the shell. Displacements are 
separately assumed for the shell, stringer, and the ring. 

TRANSVERSE SHEAR DEFORMATION FORMULATIONS 
Shell 

A consistent first order transverse shear deformation theory is developed to 
model the shell. Based on the assumption that the shell thickness t is relatively 
small and hence, does not change during loading, the displacements at an arbitrary 
material point in the shell are approximated by 


U(x,8,z) = u(x, 8) + z<j) x (x, 8) 

( 1 ) 

V(x, 8 , z) = v(x, 8) + Z<f>g(x, 8) 

( 2 ) 

W(x,8,z ) = w(x,8) 

( 3 ) 


where u(;r,0), v(x } 0) and w(x,6) are the displacements of the points of the ref- 
erence surface, and <f) x (x,6) and <f>^(x^6) are the rotations of the normal to the 
reference surface as shown in Fig. 3(a). Using Eqs. (1) to (3) and assuming small 
displacement gradients, the three-dimensional engineering strains are 


= e 7 


+ ZK 2 


Zoo 


Z$6 + ZKoq 

C 1 + it) 


e zz = 0 (4) 
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(5) 
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The transverse shear strains e ir and e# 2 represent average transverse shearing 
strains through the 'thickness of the shell since these strains contain derivatives 
of the displacements in z, and Eqs. (1) to (3) are approximate in the z-coordinat.e. 
In Eqs. (4) to (6), the two-dimensional, or shell, strain measures, which are inde- 
pendent of the z-coordinate, are defined by 


du 

^xx — 

ox 

_ d(j> x 

o 

OX 

(7) 

1 dv w 
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(8) 

dv 1 du 

1x6 = + Rdd 

(9) 

d<f>s 1 d<f> x 1 dv 

Kx9 ~ dx + R dd + R dx 

( 10 ) 

d<po 1 d(f> x 1 dv 

KxB dx R dd Rdx 

( 11 ) 

dw 

7 “ = + & 

v 1 dw 

lez = to - + £00 

( 12 ) 


If we set the (average) transverse shear strains in Eq. (6) to zero, then the 


rotations of the normal are 

dw 

Y X = 

OX 

v 1 dw 
^ 9= R~ R~dd 


(13) 

(14) 


so that 


^ x9 — ^ x9 — 


2 d 2 w 2 dv 

R dxdd R dx 


k x g = 0 


(15) 


Hence, the thickness distribution of the shear strain reduces to 
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1 x 6 + z( 1 + Yr)k x 0 
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(16) 
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which coincides with the results of Novozhilov’s (1964) classical shell theory. 

It is evident from Eq. (5) that three shell strain measures are needed to rep- 
resent the shear strain distribution through the thickness in the transverse shear 
deformation shell theory. Whereas, only two shell strain measures are required in 
classical shell theory to represent the shearing strain distribution through the thick- 
ness (refer to Eq. (16)). Also it can be shown that under rigid body rotation of the 
shell, the nine shell strain measures, given by Eqs. (7) through (12) vanish. (For 
Novozhilov’s classical shell theory, six shell strain measures given by Eqs. (7-9) and 
(15) vanish under rigid body rotations). 

The physical stress resultants and stress couples for the shell in terms of stress 
components are given, in the usual way, by 


(An, M X x) — J( 1, ^)^n( 1 d - dz 

(Noe, Mee) =J(l,z)aee dz 
(N x e, M x g) =J(l,z)a x9 (l + jj)dz 
(Ng x , M 9x ) =J( 1 , z)crg x dz 

Qx = J &xz( 1 "b jp dz 

Qe =j(rg z dz 

A generalized 9x1 stress vector for the shell is defined by 


(17) 


& shell = M X x M x q, M x $i Q X iQ$\ ( 18 ) 


in which M x q and M x q are the mathematical quantities conjugate to the modified 
twisting measures R x $ and k x q , respectively, and axe defined in terms of the physical 
stress couples by 

M x e = ~(M X 0 + Mg x ) M x g — -(M x g — Mg x ) (19) 
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The nine elements of the stress vector in Eq. (18) and the relations of Eq. (19) 
determine all the stress resultants and stress couples listed in Eq. (17) except for 
shear resultant N x e- The shear stress resultant N x $ is determined from moment 
equilibrium about the normal for an element of the shell. This so-called sixth 
equilibrium equation is 


N x $ — Nq x 


Mffx 

R 


( 20 ) 


The generalized strain vector for the shell is 

tahell = [^xx ? ^)xS' ! ^xx j KxSi Ifxz ? *)0z\ 


( 21 ) 


This strain vector is conjugate to the stress vector in the sense that the internal 
virtual work for the shell is given by 


— JJs "hell Vshell dS 


( 22 ) 


where S denotes the area of the reference surface and dS = dxRdO . This expression 
for the internal virtual work can be derived from three-dimensional elasticity theory 
by using Eqs. (4) to ( 6 ) for the thickness distributions of the strains and the 
definitions of the resultants given by Eqs. (17) and (19). 

Consistent with the transverse shear deformation theory, the linear elastic con- 
stitutive law for a laminated composite shell wall is given by 
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in which stiffnesses and Dij are given in Appendix. The transverse shear 

stiffnesses, A 44 , A 45 , and A 55 can be calculated by two different methods. The first 
method is based on the assumption of constant transverse shear strain distribu- 
tion through the thickness, and the second method is based on the assumption of 
constant transverse §hear stress distribution through the thickness. In the present 
analysis, we have used the first method to compute the transverse shear stiffnesses. 

The statement of virtual work is 


sw a n h [ u = 6w; xi + sw e x xt 


(25) 


where the external virtual work for a cylindrical shell under constant internal pres- 
sure, including an axial load due to the closed-end effect, is written as 


6W, 


ext 


SL 


p Sw dS + p 


0 

/ 

-© 


R 2 


d6 [8u(l,6)-8u(-l,d)) 


(26) 


and the external (or augmented) virtual work due to the interacting loads is 

l 

8Wx xt = J |A IS (a;)[<5u(a;,0) - -8<f> x (x, 0 )] + A zs (a;)^u)(x, 0)| dx 

-l 

& 

+ f {A rr (tf)[M0,«)-|^«(0^)]+A tfr (tf)[tfr;(0,fl)-^M0^)] (27) 

-© 

+ \„($)6w(0,e) - „„) + A„(m ^ |„ 0 

- - 9WW) - <«(-<. 0)1 

The axial force Q in Eq. (27) is an additional Lagrange multiplier that accounts 
for axial load sharing between the stringer and shell. 
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Stringer 

Stringer displacements u a (x ) and w 3 (x), and the rotation of the normal <pe s (x) 
are shown in Fig. 3(b). Based on transverse shear deformation theory, the virtual 
work expression for the stringer is 

i ' z 

j [N xa 6e xs + Me s 8 K 0 s + V Z 3 8 ~y Z3 \dx = — J |A ia (r) \ 8 u s (x) + e s 8 <pQ a (x)\ 

-i -i 

+ A Z3 (x dx + Q[dus(l) - Su a (—l)] 

(28) 

in which N xa is the axial force in the stringer, Mg a is the bending moment, V za is 
the transverse shear force, e xa is the normal strain of the centroidal line, the product 
ZK 0 S is the portion of the axial normal strain due to bending, ■y za is the transverse 
shear strain, and e s is the radial distance from the stringer centroid to the contact 
line along the shell inside surface. The strain-displacement relations and Hooke’s 


law for the stringer are 




€xS ~ M s 

K 9s — $ g 3 

Izs = 4>6s + w’ a 

(29) 

^X3 “ {■& A.') g^XS 

Mff a = ( EI) a K0 a 

V za = (GA) alzs 

(30) 


in which the prime denotes an ordinary derivative with respect to x. 

Ring 

Ring displacements are denoted u r (9), v r (d ), and w r (0), and the rotations are 
denoted by (f> xr (9 ), ), and (f> zr (0) as shown in Fig. 3(c). The structural model is 

based on transverse shear deformation theory and includes cross-sectional warping 
due to torsion. The extension of classical thin-walled, open section curved bar 
theory to laminated composite materials was developed by Woodson, Johnson, and 
Haftka (1993). However, Woodson et al. (1993) did not consider transverse shear 
deformations. 
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The statement of virtual work is 


© 

/ + M„6k xr + M zr 8K zr + T sr 8r r + M ujr 8(T r /Ro ) + V xr 8~f xr + V zr 8^ zr \* 

-© 

© 

R 0 d8 — — J | A xr {ff)\8u r (d) e r 8(f)0 r {6)\ + \e r {8) [0u r (0) + t r 8<f> xr (9) — u>o<5r r (0)] 
-© 


+ A z 


.(6)8w r (9) + Aer(0)8<p er (9) + A zr (9) [8<f> zr (9) - ^8^(8 )] } (1 + j-) Ro d8 


(31) 

in which Ng r is the circumferential force, M xr is the in-plane bending moment, 
M zr is the out-of-plane bending moment, Af wr is the circumferential bimoment, T sr 
is the St. Venant torque, V xr is transverse shear force in the x-direction, V zr is 
transverse shear force in the ^-direction, t$ r is the circumferential normal strain of 
the centroidal arc, k xt is the in-plane bending rotation gradient, k zt is the out-of- 
plane bending rotation gradient, ry is the twist rate, 7 xr is transverse shear strain 
in 9-z plane, j zr is transverse shear strain in x-z plane, e r is the distance from the 
ring reference arc to the contact line along the shell inside surface, and Rq is the 
radius of ring reference arc. Parameters u> 0 and uq are the constant coefficients in 
the contour warping function, u>(x) = o> 0 + for the attachment flange of the 
ring. The rotations and strain-displacement relations are 


1 /• 

ter = -=r(Tr + w r ) 
ft o 

T r = -ZZ-{<j>er + 4>zr) 

ft 0 
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ft o 


1 . 

7 xr — H>zr H 

Kc\ 


( 32 ) 


in which <j> xr is the rotation around x-axis, <j)g r is the rotation around 0-axis, <f> zr 
is the rotation around z-axis, and the over-dot denotes an ordinary derivative with 
respect to 8. It is assumed that the shear forces are decoupled from extension, 
bending, and torsional deformations of the ring. Thus, the material law for the ring 
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The stiffness in the first five rows and columns of this matrix were evaluated from 


a computer code developed by Woodson, Johnson, and Haftka (1993). 

For structural models in which the effect of warping of the ring cross section is 
excluded the contribution of the bimoment, M ur , to the ring virtual work in Eq. 
(31) is neglected. The fourth row and column of the stiffness matrix, Eq. (33), are 
ignored. Also, the contour warping function w(r) is taken as zero. 


CLASSICAL FORMULATIONS 


Shell 

The shell is modeled with Sanders’ (1959) theory for thin shells. Define a gen- 
eralized strain vector in terms of the shell strain measures by 

tshell = [^rn ^ 86 1 ^ 86 1 (34) 


The first five strain measures of the shell reference surface in Eq. (31) are related 
to the displacements by Eqs. (7-9), and the sixth strain measure, k x q, is given by 


d<t>e , 1 dcf> x 1 

K " - 37 + RW + R*’ 


(35) 


where the rotation about the normal, (f> z , is 

1 , dv 1 du 

4>z = Tfa-RW 

and the rotations <f> x and 4>e of the normal are given by Eqs. (13) and (14). 


(36) 
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Define a generalized stress vector in terms of the stress resultants and couples 


of Sanders’ theory by 

® shell — xxi N qq , N X Q, M xx , Mqq, M x q\ 


(37) 


such that the internal virtual work is given by Eq. ( 22 ). Quantities N* 6 and are 
the modified shear resultant and twisting moment resultant in the Sanders theory. 
Hooke’s law for a laminated composite shell wall is 


& shell = Ht s hell 


H 


A B 
B t D 


(38) 


in which the 3x3 sub-matrices A, B and D are given by classical lamination theory 
(Jones, 1975). The external virtual work expressions for the classical shell theory 
are still given by Eqs. (26) and (27), but the rotations in Eqs. (27) are given by 
Eqs. (13) and (14). 


Stringer 

The stringer is modeled with Euler-Bernoulli beam theory thereby neglecting 
the transverse shear strain. Hence, equating ~f zs in Eq. (29) to zero results in the 
following expression for (po 9 . 

<t>6s = -w' s (39) 


It may be noted that neglecting the transverse shear strain would also modify the 
virtual work statement given by Eq. (28), and the third equation in the Hooke’s 
law, Eq. (30), is neglected. 

Ring 

For classical formulations, the ring is modeled with Euler-Bernoulli beam theory 
thereby neglecting the transverse shear strains. Hence, equating 7 xr and y 2r in Eq. 
(32) to zero results in the following expressions for the rotations <f> xr and O z r . 

<t>xr — -=“(Ur ~ «>r) <j>zr = ~^~U r (40) 

ito tto 
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It may be noted that neglecting the transverse shear strains would also modify the 
virtual work statement given by Eq. (31), and the fifth and sixth equations in the 
Hooke’s law, Eq. (33), are neglected. 


DISPLACEMENT CONTINUITY 

In order to maintain continuous deformation between the inside surface of the 
shell and stiffeners along their lines of contact, the following displacement continuity 


constraints are imposed: 

Along the shell - stringer interface (i.e., —l<x<l,9 = 0), 

g xs = u{x, 0) - ^ <f> x (x , 0) - [u 3 (a:) + e s <j>e a (x )] =0 (41) 

g zs = w(x , 0) - w s (x) = 0 (42) 

Along the shell - ring interface (i.e., x = 0, — 0 < 0 < ©), 

g xr = tt(0, 9) - |<M0, 0) - [«r(«) + e r <f>8r(0)\ = 0 (43) 

ger = u(0, 9) - ^<f> 6 ( 0, 9) - [iv(0) + e r <j> xr (9) - a> o r r (0)] = 0 (44) 

g zr = u>(0, 9) - w r (9 ) = 0 (45) 

Ger = ~lix l*=° ~ = ° ( 46 ) 

G zr = [f|L=o ~ I^h=o] - lM0)-^T r (8)] = 0 (47) 

The variational form of these constraints are 

l 

f [<5A r5 g xs + 8X 23 g zs ^ dx = 0 (48) 


/[« xrQxr + fi^OrgOr “1“ 8\ zr g zr ~\~ 8AQ r Go r H - 8A. zr G zr } (R 0 + c r ) d9 = 0 (49) 
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The constraint that the elongation of the shell at 6 = 0 and the elongation of the 
stringer are the same is 

6Q{ [it(/,0) - u(—l,0)] - [u s (/) - u s (-/)] } = 0 (50) 

DISPLACEMENTS, ROTATIONS, AND INTERACTING 
LOAD APPROXIMATIONS 

The periodic portions of the displacements and rotations are represented by 
truncated Fourier Series having fundamental periods in the stringer and ring spac- 
ing. The non-periodic portions of the displacements due to axial stretching are 
represented by simple terms in x. The Fourier series reflect symmetry about the 
x-axis for the repeating unit. For the shell, displacements of the middle surface (see 
Fig. 3a) are represented as 

M N M N 

u(x,6) = ^2 u imnSin(a m x)Cos(/3 n 9)+'^2 FI u 2m nCos(a m x)Cos(f3 n 0 ) 

m = 1 n=0 m = l n = l 

(51) 

M N M N 

VimnCos(a m x)Sin(/3 n d) + EE V2 mn Sin(a m x)Sin(0 n O ) (52) 

m= 0 n= 1 m = 1 n = 1 

M N M N 

6) = EE wi mn Cos(a m x)Cos(/3 n 0) + ^ ^ w 2 m n Sin(a m x)Cos(j3 n 0) 

m = 0 n =0 m — 1 n=l 

(53) 

and rotations of the normal are 

M N M N 

4*x(.X , O') = ^ ^ ^ ^ $xlmnSlTl{oL m x)C Os( K /3 n 0^ -(- ^ ^ ^ $x2 mn^ Os[ot rn x}Cos[0 n 0 > ) 

m= 1 n=0 m—1 n=l 

(54) 

M N M N 

Ms,#) = $01 m n CoS ( Oi m X )Sin(/3 n 0) + <f>92 mn Sin( x)Sin((3 n 9) 

m=0 n = l m = l n=l 

(55) 

in which a m = and /?„ = ^ where m and n are non-negative integers. Note 
that some terms in the truncated Fourier Series of Eqs. (51-55) have been omitted. 
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The coefficients of the omitted terms are 1 / 200 , ^2m0, ^20n, 4>x 2Q0, <f>x 2 m 0 , and 

4*x20n 1 in which m = 1,2, A I and n — 1,2 , ...,7V. The rationale for their omission 
is discussed in the following sub-section. The displacements of the centroidal line 
of stringer (see Fig. 3b) are 


gif 

2 / 


uu 



M 

: + 

F u 3 imSin( 


771= 1 

M 


E 

^slm SlTl{pt m 


M 


m = l 


m = l 


M 


m=l 


and the rotation of the normal of the stringer about the 0-axis is 


(56) 

(57) 


M M 

^ ^ ^OslmSlTl^CkfnX) T ^ ^ 4 > 082mC OS^Oc m X^ (58) 

m = l m=l 

where the coeffcients u s 20 > ^520 and <f>o S 20 are omitted. Coefficient go in the axial 
displacement field of the shell and q\ in the axial displacement field of the stringer 
represent elongations of each respective element caused by either an axial mechanical 
load or due to close-end pressure vessel effects. The displacements of the reference 
circle of the ring (see Fig. 3c) are 


and rotations are 


N 

“r(«) = 5> rn Cos(fl n 6) 

n=l 

N 

v r (6) = ^2 v rnSin(ft n 6) 
n— 1 

N 

w r{Q) = w rnCos(0 n 0 ) 

71 = 0 

N 

4>xr(Q) = ^2 <t>xrnSin(P n 6) 

71 = 1 

N 

<l>6r{Q) = ^22 <t>drnCos{P n 6) 

n= 1 


(59) 

(60) 
(61) 

(62) 

( 63 ) 
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( 64 ) 


N 

<f> :r (0) = ^2 ^mSin(M) 

n= 1 

where the coeffcients u r o and <j>o r o are omitted. The distributions of the interacting 
loads, or Lagrange multipliers, are taken as 


M M 

^ ^ ><x3\mSin(0L m x) ■}- ^ ^ ( 65 ) 

rn = 1 m ~ 1 

Af M 

^ ^ d~ ^ ^ *^) (66) 

m — 1 m = 1 

N 

\ xr (9) = ^ \ xrn Cos(f3 n 9) (67) 

71=1 

N 

A<?r(0) = ^ A 0rnSin(/3 n 8) (68) 

71 = 1 

N 

Kr(6) = £ A zr n Cos(/3 n 9) (69) 

71 = 0 

N 

A er (0) = J2 A SmCos(M) (70) 

71=1 

N 

A zr {8) = A zrn Sin(fi n 8) (71) 

71=1 

where the coefficients A^sch A z2 s 0 t A xr o, and T^ r o are omitted. 


Terms Omitted in the Fourier Series 

Terms omitted in the truncated Fourier Series for the displacements, rotations, 
and the interacting loads were determined from rigid body equilibrium conditions 
for the ring and stringer, and from displacement continuity conditions between the 
shell and the stiffeners. The external virtual work for the stringer and ring must 
vanish for any possible rigid body motions of these elements. For the stringer these 
rigid body motions are spatially uniform x-direction and ^-direction displacements. 
(A rigid body rotation of the stringer in the x-z plane is not considered since this 
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motion would violate longitudinal periodicity of the repeating units.) Vanishing of 
the external virtual work for an arbitrary rigid body displacement of the stringer 
in the axial direction leads to the x-direction equilibrium equation 

/ 

J A xa (x) dx — 0 (72) 

-l 

Similarly, the equilibrium equation for a rigid body displacement of the stringer in 
the z-direction is 

l 

J A zs (a r) dx = 0 (73) 

-/ 

If the ring is considered in its entirety, that is, as made up of an integer number 
of repeating units around its circumference, the rigid body motions that lead to non- 
trivial equilibrium conditions are a displacement in the x-direction and a rotation 
about an axis through the centeroid of the ring parallel to the x-direction. The 
x-direction equilibrium equation is 

0 

J A xr (0) (Bo + e r ) d.0 = 0 (74) 

-© 

and the moment equilibrium equation about the x-axis is 

© 

J A 0r (6) ( R 0 + e r ) 2 dd = 0 (75) 

-© 

Equilibrium Eqs. (72) to (74) imply that coefficients 

Axs20 = 0 A 25 2o = 0 A xr o = 0 (76) 

in the Fourier Series for the interacting loads, and these conditions have been taken 
into account in Eqs. (65) to (67). The sine series for A$ r given in Eq. (68) satisfies 
the equilibrium condition given in Eq. (75). 
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Consider the variational form of the constraints, Eqs. (48) and (49), for the 
spatially uniform components of the virtual interacting loads. These equations are 


N 


u 200 — —<}>t200 + ^^( u 20n — n<^x20n) — (w s 20 + e s<t>Os2o) 
w n= 1 ~ 


6\ 


xs 20 


= 0 


(77) 


N 

| U>10n — 

n=0 


20 


0 


(78) 


^200 cy ^x200 “h ^ (^2m0 ^i2mo) ( ^r0 H” e r^0rO )J <^A xr o — 0 (79) 

ra=l 

Since these equations are satisfied on the basis that <5A X5 20 = 0, b\ Z s20 = 0 and 
<5A xr0 = 0, consistent with Eq. (76), the bracketed terms in Eqs. (77) to (79) do not 
necessarily vanish. The implication that these bracketed terms in Eqs. (77) to (79) 
do not vanish is that displacement continuity conditions are not satisfied pointwise. 
Pointwise continuity can be achieved by taking each Fourier coefficient appearing 
in the bracketed terms of Eqs. (77) to (79) to be individually zero. Fourier Series 
given in Eqs. (51), (54), (56), (58), (59), and (63) reflect this choice. Moreover, 
Fourier coefficients U 200 , u 5 2 o , and u r 0 represent rigid body displacement in the axial 
direction for the shell, stringer, and ring, respectively, and setting them to zero can 
be justified on the basis that rigid body displacement does not contribute to the 
deformation of the structural elements. Since Fourier coefficient w s2 0 represents 
rigid body displacement of the stringer in the 2 -direction, it would seem that it 
should be set to zero as well. However, to maintain continuity between the stringer 
and the shell in the 2 -direction, we impose the condition 


N 

E] U ’l 0 n - «>a 20 = 0 ( 80) 

n = 0 


to determine w s2 o after obtaining the solution for the displacement components 
that deform the shell; i.e., Fourier coefficients Wio n , n = 1 ,...,iV, are taken to be 
non-zero independent degrees of freedom since the stringer coefficient w s2 o is not a 
part of the solution vector. 
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Finally, consider the constraint equation associated with <f>/Wo! the spatially 
uniform component of the interacting moment intensity, which was omitted in the 
series given by Eq. (70). Derived from Eq. (49), this constraint equation is 


M 

| ^ ' a m^2ml) 


+ 


•AfJrO j 8AsrO 


0 


( 81 ) 


171 = 1 

We equated the constant component of the twist, <W o, to zero from the consider- 
ations associated with Eq. (79). Consequently, a non- zero value of the constant, 
component of the interacting moment intensity, /Wo ^ 0, would not contribute 
to the equilibrium of the ring, since /Wo and <Wo are conjugate variables in the 
external work for the ring (refer to Eq. (31)). Since <f>g r Q = 0, it is necessary that 
AorO = 0 to achieve consistent conditions for the torsional and out-of-plane bending 
equilibrium of the ring. With <!>/Wo = 0 in Eq. (81), the bracketed term does not 
necessarily vanish, and as a result pointwise rotational continuity betwen the shell 
and the ring is not assured. Pointwise rotational continuity is achieved if we take 
the coefficients u> 2 m 0 = 0, m = 1 ,...,M, as was done in the Fourier Series for the 
normal displacement of the shell given by Eq. (53). 


DISCRETE EQUATIONS AND THEIR SOLUTION 


Transverse Shear Deformation Model 

The discrete displacement vector for the shell is the (lOMiV + 3 M + 3N + 2) x 1 
vector 


- r -T ~T ~T 1 T 

V- shell — |U 0 5^1 i •••? ^A/J 


(82) 


in which subvectors are 


rp 

Uo = [?0, ^100, Viol, U^ioi, 00101? ••■? V 10N, W 10N, 0<?ion] (83) 

Um == [ w lm0 ^ilmO w 2ml?^lml^2ml)^lmli ^2 ml ? ^xlml) ^i2mli 

^lmli ^2ml) ^2mN ? 0x1 mN ? (84) 

rp> 

<f>x2mN, <j>9\mN, <f>62mN ] 
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where m — 1 , M 

The (6 M + 1) x 1 discrete displacement vector for the stringer and (6 N + 1) x 1 
vector for the ring are 


u str 



l-,U a \i,U a 2\, w aU ,w a 2i, <f>0s 21 


U'slM, W s2 M 


9slM 


, ^>0s2m] 


UslMi U s2 M, 


(85) 


T 

&ring = [^rO ? ^rl ? ^rl ? ^rl ? ftSrl ? ^rrl ? <^zrl ? •• • ? w r7V ? ^rN ? ^rN ? 4>0rN ? $xrN ? ^zriv] 

(86) 

in which the term u> s o for the stringer has been omitted as discussed in reference 
to Eq. (80). The 4M x 1 discrete interacting loads vector for the shell-stringer 
interface and (5JV + 1) x 1 vector for the shell-ring interface are 

A rji 

^ sir = [^xsll ? ^xs21 ?^zsll ? ^zs21 ? •••? ^xslM ? ^xs2M ? M ? ^zs2M ] (87) 


[A zrO i ^xrl ? ^0rl 5 ^zrl i ? -dzrl ? • • • ? ^$rN i ^zrN ? ^OrN ? -dzrTv] (88) 


Classical Model 

The discrete displacement vector for the shell is the (6 MN + 2M + 2N + 2) x 1 
vector 


- r *T - T -T 1 T 

^ shell — [Uq 5^1 ? *••? ^mJ 


(89) 


in which subvectors are 


Uo = [ 90 , W’lOO, ^101; U>101? • ••, ^ion, ^ion] 


(90) 




[^ImO j^lmOi ^lmb ^2mli ^lmb ^2 ml? ^lml ? ^2 ml 

1 T 

^2miV) ^lmN) ^2mNj 


UlmN?U2mAN 


(91) 


where m = 1 , M 

The (4 M + 1) x 1 discrete displacement vector for the stringer and (4A T + 1) x 1 
vector for the ring are 


^str 



u sll ? U S 2U w; sll ? ^s21 


W S 1M?^52M? WslM, 


W S 2M } 


(92) 
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( 93 ) 


^ ring — [^rO ? ^rl ? ^Vl j ^rl ? ^0rl i ^rjV > ^VA r ? ^r/V > ^ri\f] 

The 4M x 1 discrete interacting loads vector for the shell-stringer interface and (5 N+ 
1) x 1 vector for the shell-ring interface are the same as for the shear deformation 
model and are given by Eqs. (87) and (88). 

The approximations in Eqs. (51) through (64) for the displacements and Eqs. 
(65) through (71) for the interacting loads are substituted into the virtual work 
functionals for each structural element, and also substituted into the variational 
form of displacement continuity constraints. Then integration over the space is 
performed. (The test space of virtual displacements and the virtual interacting loads 
is the same space used for the approximations in Eqs. (5T71).) This process results 
in a lOTtf 7V + 137lf + 147V +6 system of equations for the transverse shear deformation 
model and 6 MN + 10A7 + 117V + 6 system of equations for the classical model, 
governing the displacements and the interacting loads. The governing equations 
are of the form 


"An 0 0 B\\ B 12 B 13“ &shell f An ' 

0 A 22 0 B 2 \ 0 S23 Ustr 0 

0 0 A 33 0 £32 0 ii'Ting __ 0 

Bl ; £ 2 T 1 0 0 0 0 | X str [10 ' 

B\ 2 o BJ 2 0 0 0 Xring 0 

Bl \ 0 0 0 0 J [ Q J 0 


(94) 


in which sub-matrices A'n, A’22 and A3 3 are the stiffness matrices for the shell, 
stringer, and ring, respectively. The sub-matrices Bij,i,j = 1,2,3, in Eq. (94) are 
determined from the external virtual work terms involving the interacting loads, and 
the constraint Eqs. (48) to (50). The vector on the right-hand-side of Eq. (94) is the 
force vector, determined from the external virtual work terms involving pressure. 
The constraint equations correspond to the last three rows of the partitioned matrix 
in Eq. (94). Equation (94) is first solved for the displacements in terms of interacting 
loads, then this solution is substituted into the constraint equations to determine 
the interacting loads. Thus, the total solution is obtained. 
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NUMERICAL EXAMPLE 

Data for the numerical example are representative of a large transport fuselage 
structure. The shell radius R = 122.0 in., frame spacing 21 = 22 in., and stringer 
spacing 2120 = 15.0 in. The shell wall is a 13-ply [±45, 90, 0, ±60, 90, ±60, 0, 90, ±45] T 
laminate of graphite : epoxy AS4/938 tow prepreg with total thickness of 0.0962 in. 
The ply thickness is 0.0074 in., and the lamina material properties are assumed 
to be Ex = 19.21 x 10 6 lb/in.\ E 2 = 1.36 x 10 6 / 6 /m. 2 , G 12 = G 13 = G 23 = 
0.72 x 10 6 lb/in.‘ i , and v \ 2 = 0.32. For the transverse shear deformation model, the 
shell wall stiffness sub-matrices of Eq. (23) are computed using these ply data and 
the expressions for the stiffness elements given in the Appendix. The numerical 
results are 

'0.5774 0.2619 0 

A = 0.2619 0.9766 0 x 10 6 lb/ in. 

0 0 0.2889 _ 

'3.893 0 0.1847 0.1847' 

B - 0 -5.043 -0.2213 0.2213 lb 

0 -0.443 -1.1351 1.1351 

'474.937 256.071 45.074 0 

256.071 615.194 54.003 -0.47 x 10 “ 5 

U ~ 45.074 54.003 276.965 -0.75 x 10 ~ 5 

. 0 -0.47 x 10 “ 5 -0.75 x 10 ~ 5 0.75 x 10 " 5 

and elements of the transverse shear stiffness matrix in Eq. (24) are 

A 44 = A 55 = 0.69264 x 10 5 lb/in., A 45 = 0 

The bending and stretching-bending coupling sub-matrices for classical lamination 
theory, Eq. (38), are given by 

'474.937 256.071 0 

D = 256.071 615.194 0 lb in. B = 0 

0 0 276.965 _ 

The extensional stiffness sub-matrix A is the same for classical theory and the 
transverse shear deformation theory. 
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Cross sections of the stiffeners and their dimensions are shown in Fig. 4. The 
stringer is an inverted hat section laminated from twelve plies of AS4/938 graphite- 
epoxy tow prepreg with a [±45, 02 , 90, ±15, 90, 02 , ±45] r lay up and total thickness 
of 0.0888 in. The stiffnesses in Hooke’s law for the stringer in Eq. (30) are 

(EA) S = 0.6675 x 10 7 /&, ( El), = 0.2141 x 10 7 /& in. 2 , ( GA) S = 0.843 x 10 6 /& 

We assume a 2-D braided frame consisting of 0° and 90° tows. The wall thickness 
is 0.141 inches, and the elastic modulii are assumed to be E\ = 7.76 x 10 6 lb/in. 2 , 
E 2 = 8.02 x 10 6 lb/in. 2 , G 12 = G13 = G23 = 1.99 x 10 6 lb/in. 2 , and v\2 = 0.187. 
Using the ring material properties and the cross-sectional dimensions, the stiffness 
matrix for the ring in Eq. (33) is computed. The non-zero stiffnesses are 

EA = 0.9088 x 10 7 /&, EI XX = 3.915 x 10 7 /6 in. 2 , EI ZZ = 0.1867 x 10 7 /& in 2 

EI ZX = 0.2993 x 10 7 /& in. 2 , EI^ X = -1.322 x 10 7 /& in. 3 , GJ = 0.1346 x 10 5 /fc in 2 
EI UW = 1.705 x 10 7 /& in 4 , EI UZ = -0.1865 x 10 6 /& in. 3 
GA xe = GA z e = 0.2396 x 10 7 lb 

All the results presented are for an internal pressure p = 10 psi, and the Fourier 
series are truncated at twenty-four terms in the x - and ^-directions (M — N = 24). 

RESULTS AND DISCUSSION 

Interacting Load Distributions 

The distributions of the interacting load intensities between the stiffeners and 
the shell are shown in Figs. 5 through 11. The effects of transverse shear deforma- 
tions and of warping of the ring’s cross section due to torsion on the magnitudes of 
the interacting line loads are summarized in Table 1. For the component X xs tan- 
gent to the stringer (Fig. 5), there are only small differences in the distributions as 
predicted by the four structural models. However, the peak value of the component 
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normal to the stringer, A~ 5 , is reduced in the transverse shear deformation models 
with respect to its peak value in the classical models (Fig. 6 and Table 1). 

The distributions of axial force intensity, A xr , between the ring and shell pre- 
dicted by the classical and shear deformation models with warping are nearly the 
same (Fig. 7). However, the distributions of this force intensity predicted by the 
classical and shear deformation models without warping have significant differences. 
Thus, this interacting load intensity is more sensitive to the inclusion or exclusion 
of warping of the ring cross section into the structural model. As shown in Fig. 
8, the differences in the results for circumferential force intensity, A# r , between the 
ring and shell from the four models are small, except in the vicinity of the stiffener 
intersection where the effects of including the transverse shear deformation into the 
models are manifested. However, the differences in A# r occur over one wave length 
of the highest frequency i.e., A 9/Q = 2/24. Differences occuring over the shortest 
wavelength may not be significant; more terms in the Fourier series are required to 
verify this. The distributions of the normal force intensity, A* r , between the ring 
and shell predicted by the four models are essentially the same (Fig. 9). The distri- 
butions of the circumferential moment component, Aq t , predicted by the classical 
models have higher magnitudes as compared to shear deformation models (Fig. 10 
and Table 1). Also note the change in sign of A$ v distributions in the vicinity of 
the joint as a result of inclusion of warping into the models. The classical theory 
predicts much larger magnitudes of normal moment component, A zr , compared to 
the transverse shear deformation theory for the models in which warping is included 
(Fig. 11 and Table 1). However, the reverse is true for the structural models with 
no warping. Also, there is a change in sign in the distributions of A zr for classical 
models with and without warping effects. 

The distribution of the normal component of the traction across the width of 
the attachment flange of the ring is represented by line force intensity \ zr and 
line moment intensity Aq t . The values of \ zr are nearly the same in the classi- 
cal and transverse shear deformation models (Fig. 9), but magnitudes of A# r are 
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substantially decreased in the transverse shear deformation models with respect to 
the classical models (Fig. 10). Thus, the asymmetry of the normal traction across 
the flange width of the ring is decreased in the transverse shear deformation models 
with respect to the classical models. 

The distribution of the circumferential component of the traction across the 
width of the attachment flange of the ring is represented by line force intensity A Q r 
and line moment intensity A zr . The values of \# r are nearly the same in the classical 
and transverse shear deformation models (Fig. 8). However, the magnitude of A zr 
is substantially increased in the transverse shear deformation model with respect 
to the classical model with warping excluded, and is substantially decreased in the 
transverse shear deformation model with respect to the classical model with warping 
included (Fig. 11). Thus, the asymmetry of the circumferential traction across the 
flange width of the ring is increased in the transverse shear deformation model with 
respect to the classical model without warping, and is decreased in the transverse 
shear deformation model with respect to the classical model with warping. 

For the stiffened shell configuration with asymmetrical cross section ring, the 
inclusion of transverse shear deformation and warping of ring cross section into 
the analyses influences the distributions and magnitudes of interacting line load 
components A Z5 , A xr , A$ r , and A zr . The distributions of interacting line load 
components \ xs and \ zr remain essentially the same. The cause of sensitivity to 
transverse shear deformations is two-fold: First, the tangential displacements of the 
shell along the contact lines are de-coupled from the out-of-plane rotations of the 
reference surface of the shell, and for the stiffeners the longitudinal displacements 
along the contact lines are de-coupled from the rotations of the longitudinal reference 
axes. Second, in the transverse shear deformation model, the torsional rotation 
of the ring at the shell-stringer-ring joint is de-coupled from the in-plane bending 
rotation of the stringer at the joint, thereby allowing for increased joint flexibility. 
In the classical model, the torsional rotation of the ring at the joint is constrained 
to be the same as the bending rotation of the stringer (see Fig. 12). The values 
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of these joint rotations for the four structural models are given in Table 2. Notice 
from Table 2 that the sense of the rotation changes if warping is included, and that 
the transverse shear deformation results in a torsional rotation of the ring that is 
about twice as much as the bending rotation of the stringer. 

Resultants at the Stiffener Intersection 

The interacting line load intensities acting on the inside of the shell wall can 
be resolved into a resultant at the stiffener intersection (x = 8 = 0). In general 
this resultant consists of a force with components F x , F# and F z , and a couple with 
moment components C x , C$ and C z . These components are shown in their positive 
sense on the inside of the shell wall in Fig. 13. The components of the resultant 
force vector are defined by 


/ © 


F x = 

j 

A X3 dx + / 

^ 1 
1 

dd 

(95) 


-/ 
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$ r CosO + X zr 

Sind ] ( R — 

(96) 

J 

-0 
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(R-\)d6, 


/ A 2S dx 
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/ [A zr CosO 

— Aff r Sind] 

(97) 


-/ -© 

and the components of the moment resultant of the couple are 

© 

C x = J [\ zr Sind - (1 - Cosd)X 6r ] ( R - |) 2 do (98) 

-0 

/ © 

Co — - / x A z $dx T / [ - A xr (l - Cosd)(R - |) + A er Cosd + A zr Sin6] (. R - ^ )d0 

-l -0 

(99) 

0 

C~ = J [ — A xr(R ~ -)Sind — AorSinO + A zr CosQ\ (R — — ) dO (100) 
-© 
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It is found that substituting for interacting load approximations given by Eqs. (65- 
71) into Eqs. (95-100), and performing the line integrals results in components 
F x = Fq = C x = C z — 0. Thus, at the stiffener intersection, the only non-zero 
resultants are a radial force resultant, F z , and a circumferential moment resultant 
Cq (refer Fig. 13). Iq Eq. (99) the circumferential moment component, Co, consists 
of two line integrals; first integral being the contribution of shell-stringer interacting 
loads, and second representing the contributions of shell-ring interacting loads. The 
contribution to the radial force resultant F z in Eq. (97) comes only from the shell- 
ring interacting load intensitites since the resultant from the stringer vanishes by 
Eq. (73). 

The values of the radial force and circumferential moment resultants are com- 
puted using Eqs. (97) and (99) for the four structural models under consideration, 
and are given in Table 3. The differences predicted by the four structural models in 
the magnitudes of the radial force resultant F z are very small, and are within 0.4% 
of each other. There are substantial differences in the magnitudes of circumferential 
component of the moment predicted by the four models. The values of Co predicted 
by the models with warping included are much larger than those predicted by the 
models without warping effects. The individual contributions of the stringer and 
ring to Co are also affected by the change in the model as shown in Table 3. It may 
be noted that Co is more sensitive to the effect of warping than to transverse shear 
deformation. 

Singular Behavior at the Joint 

In Table 1 the comparison of peak values of the interacting line load intensities is 
meant to convey the influences of transverse shear deformations and warping in the 
structural modeling. The peak values of components A zs , A xr , A zr , and Aq t occur 
at the joint, but these peak values do not exhibit convergence with an increasing 
number of terms retained in the Fourier series. It appears that these components 
are singular at the joint. However, the resultants F z and Co determined from these 
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line load intensities were found to converge quite rapidly. See Johnson and Rastogi 
(1994) for further discussion of this point. 

Stiffener Actions 

The distributions of the force and moment resultants in the stiffeners are shown 
in Figs. 14 through 19. The stringer axial force and bending moment distributions 
(Fig. 14) are slightly asymmetric about the origin. The bending moment distribu- 
tions in the stringer are more sensitive to the change in models as compared to the 
axial force distributions. The distribution of stringer shear force, V Z3J is shown in 
Fig. 15, and it is asymmetric about the origin. Only small differences are predicted 
by the four structural models in the distribution of V zs . 

The distributions of the circumferential force and in-plane bending moment in 
the ring are shown in Fig. 16. The differences in these distributions predicted by 
the four models are very small. The distributions in-plane shear force, V zr , in the 
ring predicted by the four structural models have negligible differences, as shown 
in Fig. 17. The out-of-plane bending moment M zr and torque T r in the ring are 
more sensitive to the change in models as shown Fig. 18. The distributions of the 
out-of-plane bending moment are symmetric about the origin, and their magnitudes 
predicted by the models with warping included are substantially larger as compared 
to the magnitudes predicted by the models without warping. The distributions of 
total torque, T r (= T sr — M ur j i?o), are antisymmetric about the origin. As shown 
in Fig. 18, the torque has reduced magnitudes in the transverse shear deformation 
model compared to the classical model when warping is included. The torque 
predicted by the models without warping is St. Venant’s torque T 5r , and this is 
negligible as shown in Fig. 18. The distributions of out-of-plane shear force, V xr , in 
the ring are shown in Fig. 19, and these distributions axe antisymmetric about the 
origin. The magnitudes of V xr predicted by the transverse shear deformation model 
are larger compared to the classical model when warping is included. However the 
reverse is true for the V xr distributions without warping. The distributions for 
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Mq s , V zs , M zt , T r and V xr shown in Figs. 14, 15, 18 and 19, respectively, indicate 
that these stiffener actions are sensitive to both transverse shear deformations and 
warping deformations. 

Shell Response 

The distribution of the normal displacement along z-curve midway between the 
stringers (0 = —0), and along the 0-curve midway between the rings (x = — /), 
are shown in Fig. 20. As depicted in this figure, there is a negligible difference be- 
tween the results from the transverse shear deformation model and classical model 
(warping of the ring is included in both models). Also, there is negligible differ- 
ence in the axial and circumferential normal strain distributions between the two 
models as shown in Figs. 21 and 22. Thus, the normal displacement and in-plane 
normal strains for the shell are not significantly influenced by the inclusion of either 
transverse shear deformations or warping deformation of the ring into the structural 
models, in part because the shell is very thin for the example studied. 

A Ring with Symmetric Cross Section 

As a benchmark for comparing transverse shear deformation model with the 
classical model, analyses were performed for a ring with symmetric cross section. 
In this case the changes made to the numerical example under discussion are to set 
the bending- coupling stiffeness EI ZX , the out-of-plane bending to warping coupling 
stiffness EI WX , and the contour warping function parameter u>o of the ring, all to 
zero. Consequently, the 0-axis, as well as the i-axis, are axes of symmetry for 
the repeating unit in terms of geometry, load, and material properties. Symmetry 
about the 0-axis implies there is no out-of-plane bending and torsion of the ring; i.e. , 
u r (0) = 4>e r (8) = 4>zr{Q) = A xr (0) = d.0 r (0) = A ZT {6) = 0 for —0 < 0 < 0. Thus, 
for the symmetrical section stiffeners only the interacting line load components 
tangent and normal to the stiffeners are non- zero. Since the internal pressure loading 
is symmetric, warping of the ring cross section does not play any role in the analyses. 
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The distributions of the tangential interacting load intensity between the shell 
and ring are shown in Fig. 23. The differences in the results from the two models are 
small except in the vicinity of the stiffener intersection. The peak magnitude of \# r 
in the transverse shear deformation model is smaller than the peak value for \g r in 
the classical model (,50.8 lb/in versus 64.5 lb/in.). However, this difference occurs 
over one wave length of the highest harmonic retained in the analysis, and may 
not be significant. The distributions of the tangential and normal interacting load 
intensities between the shell and stringer, and the normal load intensity between 
the shell and ring are not significantly different in the two models. 

For a symmetrical cross section ring, in Eqs. (95) through (100) F x = F# = 
C x — Cq — C z = 0. The only non-zero component of the force resultant is the 
radial force F z . The values of F z computed from the classical and transverse shear 
deformation models are -563.72 lb. and -561.89 lb., respectively. 

CONCLUDING REMARKS 

A unit cell model of an internally pressurized, long circular cylindrical shell 
stiffened on the inside by a regular arrangement of identical stringers and identical 
rings is analyzed. The ring is assumed to have an open asymmetrical cross section, 
and the stringer is assumed to have a symmetric section. The asymmetrical section 
ring significantly complicated the analysis of the unit cell, since symmetry about 
the plane of the ring is lost. Out-of-plane bending and torsion of the ring occur 
as well as a rotation of the shell-stringer-ring joint about the circumferential axis 
of the ring. Mathematical formulations for the linear elastic response presented in 
this work include the effects of transverse shear deformations and of out-of-plane 
warping of the ring’s asymmetrical cross section due to torsion. Closed-end pressure 
vessel effects are included in the analyses. Data representative of a large transport 
aircraft are used in the numerical example. 

For the stiffened shell configuration with an asymmetrical cross section ring, the 
inclusion of transverse shear deformation and warping of ring’s cross section into 
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the analyses influences the distributions and magnitudes of interacting line load 
components between the shell wall and stiffeners (Table 1), and the stiffener actions 
(see Figs. 14 and 18). However, the normal displacement and in-plane normal 
strains for the shell are not significantly different in the four structural models. The 
cause of sensitivity to transverse shear deformations can be attributed to the de- 
coupling of the torsional rotation of the ring at the shell-stringer-ring joint from the 
in-plane bending rotation of the stringer at the joint, thereby allowing for increased 
joint flexibility. In the classical model, the torsional rotation of the ring at the joint 
is constrained to be the same as the bending rotation of the stringer (see Fig. 12 
and Table 2). The inclusion of warping of the ring’s cross section due to torsion into 
the analyses causes the magnitude of the circumferential component of the moment 
resultant of the interacting line loads, Co, to increase substantially (see Table 3). 
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Appendix 

ELEMENTS OF STIFFNESS MATRIX FOR A CYLINDRICAL SHELL 
BASED ON TRANSVERSE SHEAR DEFORMATION THEORY 

Based on the transverse shear deformation theory, the elements A t j, B t j, and 
Dij of the stiffness matrices, in Eqs. (23) and (24) for the constitutive law for a 
laminated shell wall, are given by 
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where Qij are the transformed reduced stiffnesses given in the text by Jones (1975). 

Based on the assumption of constant transverse shear strain distribution through 
the thickness, the transverse shear stiffnesses are given by 
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C 45 =(Gi3 — G2z)CosaSina 

C55 —G23C0S 2 a + Gi^Sin 2 a 
in which a is the ply orientation angle. 

Based on the assumption of constant transverse shear stress distribution through 
the thickness, the transverse shear stiffnesses are given by 
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Table 1: Effect of transverse shear and warping on the interacting load intensities and their 

distributions along the contact lines 8 



Peak values of the interacting load intensity 


Component 

CLT b 

model 

Warping c 

neglected 

CLT model 
Warping 
included 

SDT d 
. model 
Warping 
neglected 

SDT Model 
Warping 
included 

Comments on the 
distribution, and 
figure number 

^xs 

lb/in. 

131. @ 
x/l = - 0.22 

132. @ 
x/l = - 0.22 

127. @ 
x/l = - 0.21 

130. @ 
x/l = - 0.22 

Antisymmetric; 
Non-zero over 
entire stringer 
length; Fig. 5 

^ZS 

lb/in. 

767. @ 
x/l = - 0.02 

613. @ 
x/l = 0.01 

571. @ 
x/l = 0 

571. @ 
x/ 7 = 0 

Symmetric; Small 
magnitudes except 
near origin; Fig. 6 

^xr 

lb/in. 

90.8 @ 
6/0 = 0 

74.6 @ 
0/0 = 0 

10.7 @ 
0/0 = 0 

80.7 @ 
0/0 = 0 

Symmetric; Small 
magnitudes except 
near origin; Fig. 7 

lb/in. 

63.2 @ 

0/0 = ±0.04 

55.2 @ 

0/0 = ±0.04 

49.0 @ 

0/0 = ±0.21 

44.8 @ 

0/0 = ±0.29 

Antisymmetric; 
Non-zero over 
entire ring length; 
Fig. 8 

^zr 

lb/in. 

-886. @ 
0/0 = 0 

- 883. @ 
0/0 = 0 

- 854. @ 
0/0 = 0 

- 852. @ 
0/0 = 0 

Symmetric; Small 
values except near 
origin; Fig. 9 

A 0r 

in.-lb/in. 

- 198. @ 
0/0 = 0 

73.1 @ 
0/0 = 0 

- 30.2 @ 
0/0 = 0 

22.1 @ 
0/0 = 0 

Symmetric; Nearly 
zero except near 
origin; Fig. 10 

A zr 

in.-lb/in. 

3.21 @ 

0/0 = ±0.03 

-4.38 @ 
0/0 = ±0.37 

- 3.95 @ 
0/0 = ±0.37 

- 1.55 @ 
0/0 = ±0.20 

Antisymmetric; 
Non-zero over 
entire stringer 
length; Fig. 1 1 


a. Results for Fourier series truncated at 24 terms in the axial and circumferential directions. 

b. CLT is classical lamination theory. 

c. Out of plane warping of the ring’s cross section due to torsion 

d. First order transverse shear deformation theory 



Table 2: Rotations about the circumferential axis at the stiffener intersection. 


Description of the rotation 
of the structural component 


Rotations in 10' 5 radians 


Classical theory 

Transverse shear theory 

No Warping 

Warping 

No Warping 

Warping 

Shell normal ^ (0, 0) 

-2.56 

2.58 

-1.06 

2.65 

Ring twist <t> 0r (0) 

-2.56 

2.58 

-2.67 

3.64 

Stringer normal <t> ef (0) 

-2.56 

2.58 

-0.29 

1.85 


Table 3: Resultants at stiffener intersection. 


Components of 
the resultant 

Classical theory 

Transverse shear 
deformation theory 

Warping 

No warping 

Warping 

No warping 

C e from 
stringer, lb-in. 

-0.0921 

- 1.1696 

- 0.7797 

- 0.2953 

Cq from ring, 
lb-in. 

5.645 

1.627 

6.0192 

1.363 

Cq total, lb-in. 

5.5526 

0.457 

5.2396 

1.0676 

F z , lb. 

- 564.56 

- 564.06 

-563.15 

- 562.27 








(a) 




Fig. 3. Displacements and rotations for (a) shell, (b) stringer, and (c) ring 
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Stringer 

a = 1.4 in. 
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Ring 

h = 5.6 in. 
b = 1.4 in. 
e r - 2.854 in. 


Fig. 4. Stiffener cross sections. 
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Fig. 5 Stringer-shell tangential force intensity 
in axial direction. 


(N/mm) 



(Ib/in) 


<< 


w 

N 


Classical theory 

Shear deform, theory 

Classical theory w/o warping 



140.1 

105.1 

70.1 
35.0 
0.0 


- 35.0 


- 70.1 


-1 - 0.5 0 0.5 1 


x/l 



Fig. 6 Stringer-shell normal force intensity. 
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Fig. 7 Ring-shell tangential force intensity in axial direction. 
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Fig. 8 Ring-shell tangential force intensity 
in circumferential direction. 
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Fig. 9 Ring-shell normal force intensity. 
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Fig. 10 Ring-shell tangential moment intensity 
in circumferential direction. 
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Fig. 11 Ring-shell normal moment intensity. 
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Fig. 12 Joint rotations for (a) the classical structural models, and 
(b) for transverse shear deformable models. 
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Fig. 14 Stringer axial force and bending moment 
distribution. 
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Fig. 15 Stringer transverse shear force distribution. 
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Fig. 16 Ring circumferential force and in-plane 
bending moment distribution. 
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Fig. 18 Ring out-of-plane bending moment and torque distributions. 
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